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Abstract 

This paper aims to explore the inherent connection among Heisenberg 
groups, quantum Fourier transform and (quasiprobabihty) distribution 
functions. Distribution functions for continuous and finite quantum sys- 
tems are examined from three perspectives and all of them lead to Weyl- 
Gabor-Heisenberg groups. The quantum Fourier transform appears as 
the intertwining operator of two equivalent representation arising out of 
an automorphism of the group. Distribution functions correspond to cer- 
tain distinguished sets in the group algebra. The marginal properties of a 
particular class of distribution functions (Wigner distributions) arise from 
a class of automorphisms of the group algebra of the Heisenberg group. 
We then study the reconstruction of Wigner function from the marginal 
distributions via inverse Radon transform giving explicit formulas. We 
consider applications of our approach to quantum information processing 
and quantum process tomography. 

1 Introduction 

Quasiprobability distribution functions (or simply distribution functions) on a 
quantum system provide an alternative and equivalent description of quantum 
states. We will discuss three possible approaches to distribution functions. The 
first approach is essentially Wigner 's original approach [T] and it attempts to 
give a "phase-space" description of quantum states. The state of a quantum 
system determines the probability distributions of its observables. It is possi- 
ble to completely specify the state by giving the distributions of functions of 
a pair of conjugate nondegenerate observables / and g. The most well-known 
example is the position-momentum pair. Thus corresponding to a (mixed) state 
p we associate a real function W(j), q : p) oi "c-number" variables which have 
the same information content as the state and give the correct marginals. Now 
for conjugate observables the expectation values of the functions 0(/,g) will 
specify the state completely (ignoring the questions of operator ordering). In 
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classical probability theory these expectation values are generated by the char- 
acteristic function. The characteristic function of (classical) random variables 
Xi,X2, . . . , Xn with joint probability distribution . . . , Xn) = F{-x) is given 

by m 

F{t) = Je**''dJ^(x) 

where t-x is the usual Euclidean scalar product of two real vector x = (xi , • • • , Xn)^ 
and t = (ii, • • • where T denotes transpose. If the probability distribu- 

tion is given by a probability density J5(x) then the characteristic function is 
simply the Fourier transform of p(x). Another way of viewing the characteristic 
function is to note that F{t) = (e'*'-'^), the expectation value of the complex 
random variable e** -^. Then assuming the existence of a probability density it 
is given by the inverse Fourier transformation 

p(x) = Je"'*"<e**-^>dt (1) 

In the form ([T]) it is suitable for a "quantum" extension |3]. Of course, in the 
quantum case the observablcs (or quantum random variables) Xi will not com- 
mute in general and we have the problem of interpreting the function p as a 
joint probability distribution. However for a set of compatible or commuting 
observablcs a joint distribution is unambiguously defined. For incompatible ob- 
servablcs wc may take ([1]) as the definition of joint probability distribution. The 
values that are obtained on joint measurement of these observablcs constitute 
the joint spectrum which, in general, may have both continuous and discrete 
segments. In the finite-dimensional case we have a finite spectrum and hence 
the integral has to be replaced by a sum. But there are problems in inter- 
preting this as a function in a classical phase space Alternatively, we can 
work with finite Fourier transform. Such a transform is defined over a finite 
abelian group. From the structure of such groups we may focus on the group 
Zjq = TLjXI, the additive group of integers modulo iV. Wc can thus take our 
"configuration space" to be Zjv- This in turn forces us to consider observablcs 
which take "values" in the set {0,1,...,A^ — 1} where N must now be identi- 
fied with the dimension of the Hilbert space. Since we are concerned only with 
probability distributions of the possible outcomes this is really not a restriction. 
Operationally, we can always calibrate our instruments to yield these outcomes. 
The "position" and "momentum" variables are both identified with Zjv and the 
corresponding unitary representations respectively act multiplicatively and ad- 
ditively on the "position space" . So we have two representation of Zjsi. But they 
do not constitute a representation of the "phase space" Zj^ x Z^r as the latter 
is a commutative group. The simplest possible noncommutative extension is a 
central extension [5] . After some restrictions due to finitcncss of the dimension 
we get a Heiscnbcrg group. 

Now let us approach the problem of distribution function from another per- 
spective. The state of the system, p is a positive definite operator with trace 1. 
In infinite dimensions they belong to a special class called trace class operators. 
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In finite dimensions every operator is clearly trace class. Trace class opera- 
tors admit a Hilbert space structure. Thus for two such operators A, B define 
{A,B) = Tr(A^i?). In finite dimensions this introduces the familiar Frobenius 
or Hilbert-Schmidt norm. If we restrict to hermitian operators we get a real 
Hilbert space K.. Pick an orthonormal basis {Ai : i = 1, 2, . . . , } in A^. We can 
write the state 

P = 2 ^ith W,{p) = Tr(pA,) 

i 

Clearly Wi are real. If we also demand that Wi = 1 then we have a quasiprob- 
ability distribution over the index set I {i e I). We want this index set to have a 
classical interpretation and a natural choice is the phase space. Then, i = (x, z) 
is a pair of indiccfl Henceforth, wc assume this and write A(x^ z) instead of i. 
The works [SI E] follow this approach to distribution function (sec also jH| for a 
review in the finite-dimensional case). Thus the choice of distribution is equiv- 
alent to choice of special bases. The operators A{x, z) arc called phase-point 
operators. Actually, they have to satisfy some extra conditions. Wc will see 
that the phase point operators correspond to certain sets (called Wigner sets) 
in the group algebra of the Heisenberg group. 

There is yet another view of distribution function which has origins in signal 
analysis. A signal may be represented in the time domain as /(<) or frequency 
domain as /(w) (Fourier transform). Thus we represent the signal in terms of 
"elementary" harmonic signals and the coefficients give the representation in 
frequency domain. But it can also be represented by other elementary non- 
harmonic signals with minimum uncertainty. This was Gabor's seminal idea 
[S]. Unlike harmonic signals Gabor's elementary signals arc localized in time 
and frequency domains. This joint time-frequency domain is the analogue of 
phase space. How do wc generate these elementary signals? Starting with a 
"reasonable" initial signal say a Gaussian function in the time domain we apply 
a sequence of two operators, multiplication and translation, which are "diago- 
nal" in the time and frequency domain respectively. The resulting sequence of 
functions are used to represent an arbitrary signal. Now to represent a vector 
in some space we don't need a basis, any set of vectors that span the space will 
do. In finite dimensional Hilbert space such sets define frames jlO| . An exam- 
ple of such overcomplete sets is the set of coherent states in quantum optics. 
The frame-theoretic approach to distribution functions was recently proposed 
in [11) . We will not go into details of frame theory but mention that frames are 
a generalization of orthonormal bases in Hilbert space. In this context, Gabor's 
elementary functions constitute Gabor-Weyl-Heisenberg (GWH) frames. Most 
distribution functions are examples of such frames although there are some ex- 
ceptions. The GWH-frames are generated by applying sequence of translation 
and multiplication operators to (continuous) signals creating a function in the 
time-frequency domain (the phase space!). These operators generate a discrete 
Heisenberg group. 

^ We often use z in place of p for momentum variable to avoid confusion with probability 
density. Further, x and z can be vectors. 
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Finally, we come to another significant property of the Wigner function, 
a particularly important distribution function. Let W{x,p) be the continuous 
Wigner function with x and p representing the classical position and momentum 
variables. Then the marginals W{x,p) and W{x,p) are probability distri- 
butions of the quantum momentum and position operator respectively. Further, 
if we sum W{x,p) along some line ax + bp = then the resulting function is the 
probability distribution of a quantum operator "orthogonal" to cx + dp where 
(c, d) is a vector orthogonal to (a, 6) in the x — p plane. We will make these 
definitions precise later. Call this the Radon property. This is an important 
property and can be used to invert the transform. We will see that the Radon 
property is related to the transformation properties of distribution function un- 
der some automorphisms of the Heisenberg group. We note that the Radon 
property is very useful in the practical problem of reconstruction of states and 
processes. We prove a general Radon property which gives us a lot of freedom 
in our choice of possible measurements. 

The brief (and incomplete) survey of the approaches to distribution functions 
in the preceding paragraphs indicates that a lot of work has been done in this 
aregd. Besides their theoretical significance distribution functions have appli- 
cation in state tomography [T3], statistical mechanics and quantum optics [13]. 
It is also intimately connected with the theory of coherent states. The GBH 
type operators after complexification and some algebra give rise to the familiar 
displacement operators. The coherent states are the orbits of Weyl-Heisenberg 
group (henceforth only Heisenberg group) [15] . In this work we mainly focus 
on the finite-dimensional case. This case presents some difficulties absent in 
the continuous case. The finite dimensional case is also significant for quantum 
information processing [16l [17] . 

We have considered Zjv x ^at as the basic model of finite "phase space" . In 
the literature other phase spaces have been considered (see [11] for a discussion 
and references). The intrinsic structure of these phase spaces may have inter- 
esting bearing on the corresponding distribution functions. In particular, some 
authors have considered finite field F]\[ with N elements instead of (see e.g. 
[T5]). This is only possible if = is a power of some prime p. Now Zp 
and Fp coincide. In general, the additive groups of FpTz and Zp x ■ ■ ■ x Zp {n 
factors) are isomorphic. We can consider the Heisenberg groups over the latter 
(by central extension). We do not follow this here as the paper is already quite 
large. However, note that the analogy between what the authors in [18] call 
quantum nets and Wigner sets defined in this paper. More precisely, quantum 
nets correspond those Wigner sets which are permuted by the action of the au- 
tomorphism group SL{2, Zn)- The paper is quite self-contained. We give most 
of the proofs. Some of the results are known but were derived using different 
approaches. Let us first note some of the main contributions of the present 
work. 

1. We use the Heisenberg groups as the basic approach to distribution func- 
tions. As we have seen in the preceding paragraphs this is the unify- 

^See |12| and [S] for large lists of references relevant to the present work. 
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ing thread tying the different approaches and perspectives on distribution 
function. 

2. Hcisenbcrg group has been used in the hterature context of distribution 
functions. But here we use discrete Heisenberg groups and family of finite 
quotient groups thereof. We define these groups abstractly in terms of 
generators and relations. Thus we can consider different representations 
(irreducible and reducible) and operators between representations. 

3. Our treatment of the Heisenberg groups is defined in terms of generators 
and relations. This makes the computations and proofs easier. Further, 
we don't need the language of (pseudo) phase spaces. 

4. We show that the existence of distribution functions is equivalent to cer- 
tain sets in the group algebras of the Heisenberg groups. This provides 
us with powerful methods of representation theory. We list some of the 
outcomes by the use of these methods. 

(a) The analysis of marginal distributions become transparent. They 
correspond to an invariance up to permutations under certain groups 
of automorphisms (e.g. SL(2, Zn)) of the Heisenberg groups. 

(b) Demanding this invariance we get unique distribution functions in 
odd dimensions. 

(c) We also infer that it is impossible to retain full invariance and other 
properties like hermiticity and linear independence even dimension. 
We therefore have three possible strategies: i. drop the requirement 
of invariance, ii. drop the requirement of independence or hermiticity 
or iii. require invariance under a smaller set of transformations. We 
discuss all three and give some alternative candidates for distribution 
functions in even dimensions. 

(d) Our analysis via the automorphism groups obviates the need for ad 
hoc hypothesis and guess work. 

We note again that although the results mentioned in (a), (b) and parts 
of (c) are known our approach via automorphisms is different. 

5. We give explicit formulas in most cases. The close connection with finite 
Radon transform is made clear. It is used to derive the formulas for state 
reconstruction. The inversion formulas in the case o dimension = 2" the 
formulas appear to be new. 

6. We explore applications to quantum computing and information. The fact 
that Weyl-Heisenberg groups describe the kinematics of quantum systems 
is known jl9[ ?IU\ . We show that the dynamics is described by (unitary) 
automorphisms of the group algebra. The case of unitary automorphisms 
of the group itself is analyzed in [3T] . The latter correspond to the Clifford 
group and to go beyond it ( "non-classical" dynamics) we have to consider 
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the group algebra. We also illustrate the utility of the Heisenberg groups in 
quantum circuits. We give an application to quantum process tomography. 

We now give a summary of the paper. In Section [2] wc discuss the quantum 
Fourier operators. The quantum Fourier transform (QFT) is another form of 
finite Fourier transform |22j . Consider the two representations of acting as 
multiplication and translation respectively. The QFT connects the two. The 
generators of the two representations give us the basic operators: Z and X. 
We can consider these as the finite-dimensional analogue of unitary operator 
generated by "position" and "momentum" respectively In Section [3] we re- 
view (continuous) quasi-probability distribution functions or simply distribution 
functions. The continuous distribution functions are somewhat easier to deal 
with because the "infinitesimal" generators satisfy simple commutation rules 
(the Heisenberg relations). In Section |4] we come to one of our main themes the 
finite distribution functions. We list a set of properties, satisfied by the contin- 
uous Wigner function, and demand that any distribution function must satisfy 
them. In particular, we give examples of discrete Wigner functions. Here we 
encounter the difficulties when the dimension is even. We also derive explicit 
formulas for the phase-point operators. The odd-dimensional case (apart from 
some constants) is essentially same as Wootters' |7j operators in prime dimen- 
sion. 

In Section [5] we study the Heisenberg groups. We start with the continuous 
version as it has been studied well in connection with Fourier transforms |23| . 
We then look at discrete and finite Heisenberg groups, their structure, repre- 
sentation and automorphisms, all of which play important role in our study of 
the finite distribution functions. We show that there is a one-to-one correspon- 
dence between distribution functions in dimension N and certain sets {A{x, z)}, 
(called Wigner sets) in the group algebra of the Heisenberg group H in that 
dimension. The representation of these sets are the phase-point operators. A 
slight generalization of the Wigner sets may be used to define Weyl- Heisenberg 
frames. The group SL{2, Zm) of 2 x 2 matrices in Zm with determinant 1 induce 
automorphisms on the Heisenberg group H. Thus for each M G SL{2, Zm) we 
define an automorphism of H. These automorphisms, in turn, determine 
the marginal properties of the Wigner functions. Thus if W{x, z) = W{C,) is a 
distribution function then the functions 

Q{z) =Y,W{M~\) and P{x) =Y,W{M-^C) 

X z 

are the marginals. Q{z) is called a simple marginals if it is the probability 
distribution (in the given state) of an observable zjv/ defined by e*^*^ = gm^Z). 
A similar definition can be given for P{x). We show the necessary and sufficient 
condition for the existence of simple marginals for all members of SL{2, Zf^). 
Thus the requirement that all marginals be simple determine the distribution 
function up to isomorphism. In the case of odd dimensions for the Wigner 
function all marginals are simple. This can be neatly expressed as follows. 
Let Am{x,z) = A{M~^x, M"^ z). Then {Am{x, z)} is also a Wigner set. An 
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analogous result is called Clifford invariance in This is not true in even 

dimensions. We investigate three alternatives by weakening our requirements. 
First, we do not demand that the phase-point operators form a basis. Now 
they constitute a frame. This is the most common approach (see for example, 
PSI [TB]). The marginal conditions are simple but at the expense of losing 
orthogonality of bases. We show that this is similar to the case of spin-1/2 
representation in the sense that a complete "rotation" does not preserve the 
values of the functions involved. More precisely, we get functions which are not 
periodic on Z^- However, they have period 2N. So we go over to Z2N x Z2N 
as the phase space. Next we drop the requirement that the marginals be simple 
in the above sense. It is still possible to compute the marginals in terms of the 
probability distribution of the observable zm- We indicate how explicit formulas 
can derived to compute this. Finally, since it is not possible to satisfy simple 
marginal conditions on all of 5^(2, Zn) we consider certain subsets adequate for 
inversion, that is, computing the state from the marginal data. We give such a 
subset in dimension N = 2^ . A similar construction from a different perspective 
was done in [26| but we present our formulas in an explicit functional form. 

In the section [01 we explore the fact that the definition of marginals is a 
Radon transform of W in the sense of [571 HSl HH] ■ We then give the inversion 
formulas in several cases. The inversion formulas for odd dimensions were given 
in |13] . Our derivation, however, is more general and applicable to any finite 
distribution function. The important point is we can invert these transform 
and recover the Wigner function and hence the quantum state. In the next 
section we discuss some applications to quantum information processing. We 
provide some simple relations between standard quantum gates and operators 
representing the Heisenberg groups which will prove useful for implementing the 
state and process determination schemes using phase-point operators. We give 
formulas for quantum process tomography using phase-point operators. In the 
final section wc discuss some more possible applications and future work. 

2 Quantum Fourier Operators 

Let G = Z/NZ be the additive group of integers modulo- A^. There arc two 
obvious representations of G on an A^-dimensional Hilbert space H. Let g be a 
generator of G. Suppose (j) : G ^ U{H) is faithful representation of G by unitary 
operators where U{H) is the set of unitary operators on H. If 4>{g) = Z then we 
must have Z^ = 1 since the order of G is A' and the representation is faithful. 
The eigenvalues of Z are A'th roots of unity. Let {|z):i = 0,---,A^— l}be the 
corresponding eigenvectors such that Z\i) = |«) where w is a primitive A^th 
root. Call it the the computational basis Be- There is another representation 4/ 
of G defined by (t>'{g) = X where X\i}=\i+l (mod A^)>. We can think of </> 
as the multiplicative and 4>' as the additive representations. Z and X represent 
multiplication and translation operators resp. Clearly, X is unitary and there 
is basis in which it is diagonal. The quantum Fourier transform (QFT) 
is the unitary map connecting the two representations taking Be — > Bf. The 
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eigenvalues of X are also roots of unity as = I. Since is also faithful the 
diagonalization of X yields Z fixing the ordering. Hence there exists a unitary 
operator fl such that 

n^'xn = Z (2) 

The explicit form of fl in the computational basis is easy to compute. Thus, if 
a = J]^Xi \i} is and eigenvector then Xa = ua implies xq = uxi, . . . ,xn-2 = 
uxn-1, and xn-i = uxq. This yields after normalization 

m3 = -4^^-'' (3) 



So the quantum Fourier transform (QFT) is the map 



3 3 

We note that we follow the convention of mathematicians in the definition of 
discrete or finite Fourier transform. In the quantum information literature the 
usual definition is with a positive sign in the exponent which is our inverse 
Fourier transform. Now X and Z can be expressed as 

X = and Z = e" (4) 

where x and z are the hermitian generators of the respective unitary rotation. 
Moreover, their eigenstates are (discrete) Fourier transforms of each other. This 
is reminiscent of position and momentum observables which also have the prop- 
erty that their (generalized) eigenstates are (continuous) Fourier transforms of 
each other. We may therefore regard the observables x and z as conjugate 
"dynamical variables" . This terminology is further justified by the following 
observation which is crucial for our calculations of quasi-probability distribu- 
tions. 

XZ = ujZX (5) 

This is most easily derived by applying both sides to vectors in the compu- 
tational basis Be- We observe that the unitary operators e^°-P and e'^'' corre- 
sponding to translations in (continuous) position and and momentum space 
respectively obey a similar relation. 

Suppose now that H is a, product space, that is, H = ®"^Hd where Hd is a 
d-dimensional space and N = d™-. As a simple application of the basic relation 
([2|) we show that the Fourier transform of product states in the computational 
basis arc also product states and generalize a computationally useful formula. 

Lemma 1. // the basis Be, the eigenvectors of Z , consists of m-fold product 
states then their Fourier transforms are also product states given by 

m— 1 m—1 m—1 

nm = ( \r}) ® • • • ® ( 2 oj-^'"- |r» ® ( co-^^ |r» 

r=0 r=0 r=0 
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Proof. Observe that there is an impUcit ordering of the product states. Thus if 
J = 1jT=o ^^^r is the representation of a positive integer ^ j ;g d™ — 1 then 
the state \j) = |jm-i> (g) • • • (g) |ji> ® | jo> = |jm-i> • • • I ji> |.7o> = | jm-i • • ■ io> 
where we suppress the tensor product symbol in the last two relations. Further, 
we write |0> for |0---0> 

From the definition of QFT 

(km-l ■ ■ ■ fco|f^|j,„_l • ■■jo) = (0\X-''n\j,n-i ■ ■■jo) 

= (o\nn^x~''n\j^^i ■ • • jo> = (o\nz~''\jm-i ■ ■■jo) 



'<o|r!|j,„_i---jo> 



A direct computation shows that for the state 

^ m — 1 m — 1 m — 1 

|^,> = ^( 2 LO-^'"^''^ |r» ® • • • ® ( 2 \r)) ® ( 2 c.-^'- |r» 

<A:|V'j> = w-'^J. Since this is true for aU ^ A: iV - 1, n\j) = \ipj). □ 



3 Distribution functions in quantum systems 

One of the motivating factor's for distribution functions in Wigncr's work [1] 
was the construction of a quantum analogue of Liouville density in classical 
phase space. Following this approach suppose we want a "joint" distribution 
function of the operators X and Z. More precisely, we seek hermitian operators 
X and and z such that 

X = and Z = e'^ (6) 

and then try to find distribution functions associated with the observables x 
and z. We will do our computations in the computational basis Be in which Z 
is diagonal. It is easy to find z. Thus 



2tt 







Vo 



(7) 



iV- 1/ 



Of course z is only determined modulo 27rfc. From ^ and ^ we have 



{x)ij = (r2^zr2)ij = ^ kuj 



Tr{N - 1) if i = j 

2tt N-1 



k=0 



N 



(8) 



These entries of the matrix x imply that it is a hermitian circulant matrix. But 
a general linear combination ux + vz is not a circulant but a Toeplitz matrix. 
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There are efficient algorithms for finding the eigenvalues and eigenvectors of 
such matrices. So in principle we can compute expressions like ^e'^"^"*""^^) for 
real or integer u, v using the standard diagonalization procedure. It is feasible 
to find analytic expressions, however, in low dimensions only. We will tackle the 
problem by different approach. Let us briefly review the continuous case first. 

3.1 The Wigner distribution 

In the case of canonically conjugate variables like position and momentum a 
number of quasi-probability distributions are possible, each corresponding to 
a an operator ordering prescription. This is facilitated by the fundamental 
commutation relation 

[q,p] = ih 

between the position and momentum operators. Taking traces it is clear that 
such a relation is not possible in finite dimensions. So in finite dimension it is not 
clear how to prescribe ordering of operators. Moreover, there is some ambiguity 
in defining hermitian generators themselves. For example, for integers a and 6, 
x' = X + 2'KaI and z' = z + 27r&/ are also infinitesimal generators for X and 
Z respectively but their linear combinations give rise to different set of unitary 
operators. The problem of non-uniqueness is essentially the same as the one 
that arises in defining roots and logarithms of complex numbers. Hence, we 
restrict to the principal branch of the logarithm as evident in the definition of 
X and z. 

3.1.1 Continuous Wigner distribution 

The inversion formula of a characteristic function of classical probability is differ- 
ent for continuous and discrete probability distributions. In finite-dimensional 
quantum systems or more generally in the discrete part of the spectrum of a 
quantum observable we should use a formula analogous to that for discrete dis- 
tributions [5]. But, this gives a quasi-probability distribution which may not 
have the desired properties [4|. The problem seems to be rooted in the non- 
commutativity of quantum observables. The continuous Wigner distribution is 
defined by 

^) = J J<e*("*+"^)>e-'(""+"")dudv (9) 

In this equation and the rest of the paper, unless the limits arc explicitly stated, 
the real integrals arc over the whole real line. Further we use the notation 
r = (x, y)^ for a real vector in 2 dimensions. The following theorem gives some 
of the important properties of the continuous Wigner distribution. First, we 
make the dependence on the state p (mixed state, in general) explicit when 
necessary: Wc{x, z : p). We give a simple proof of a well-known results in the 
appendix. 
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Theorem 1. The function Wc{x,z : p) is real. Moreover, it gives the correct 
marginal distributions. 

dz \Wc{x,z : p)dx = \ (z\p\z}dz (10) 

where jz) are generalized eigenvectors of z. We have a similar relation for 
the other marginal. We also have the following results on general marginal 
distribution. Let R be an orthogonal matrix of order 2 representing a rotation. 
Let 

= Rr with R = ( '"'^^ „ ^™ 
sint^ cos( 

Similarly we define orthogonal transformation in the "noncommutative" space. 
If 



x'\ _ j^fx\ _ / cosOx + sinOz 
z' / \z } \ — sm9x + cos 9z 



The 



dz' \Wc{x,z: p)dx' = \ dz'{z'\p\z') (11) 

J a J J a 

where |z') are generalized eigenvectors of z' and the variables x and z are con- 
sidered as functions of x\ z' . A similar result holds for the conjugate observable 

X. 

The continuous version of distribution function of discrete observables is 
problematic. First, wc want the marginals to resemble classical discrete proba- 
bility distributions so we have delta functions at the isolated points. To justify 
the later we have to integrate over some domain of the continuous variable 
and this causes some problems interpreting these as probability distributions. 
Some authors have attempted to tackle these problems by focusing on continu- 
ous families of discrete observables like spin direction. These approaches seem 
somewhat unnatural to us. Discrete distributions must be characterized by a 
discrete measure (e.g. the counting measure) and thus the integrals must be 
replaced by sums. In particular, for finite systems we must have finite sums. 
This is the avenue we will explore in this paper. Finally, let us mention an 
important point. The Wigncr distribution and some other related probability 
and quasi-probability distributions are sometimes interpreted as joint probabil- 
ity distributions of incompatible observables. Clearly, any measurement of such 
distribution must give unsharp values of these observables, consistent with the 
uncertainty principle. The Arthrus-Kelley scheme [301 [31] is an example. For 
discrete observables the concept of joint distribution of noncommuting observ- 
ables is difficult even for fuzzy measurements. 



4 Discrete Quasiprobability Distributions 

The Wigner and other distribution functions are an alternative to the density 
matrix formulation of quantum theory and are given by distribution function 
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W{y : p) with y representing classical parameters. Expectation values of any 
quantum mechanical quantity that can be computed in a given state p can 
be computed from W(y : p). Hence we have a correspondence between quan- 
tum observables and "classical" observables along with an ordering prescription. 
Since, the density matrix provides a maximal description of a quantum system 
so does W{y : p). We thus have an alternative scmiclassical picture. In some 
situations the latter may be easier to determine experimentally In any case, 
such quasiprobability distributions provide a useful tool for visualization. 

4.1 Properties of Distribution Functions 

Let us make precise the requirements wc impose on distribution functions. Let 
W{y : p) be a distribution function associated with a quantum state p of a 
quantum system S and y is a real vector representing a finite set of "classical" 
parameters. Let H be the system Hilbcrt space and S{H) the set of states, that 
is, the convex set of normalized positive trace-class operators. 

Rl. W{y : p) is a continuous real function on S{H) that preserves convex 
combinations: if pi, p2 e S{H) and s ^ 1 then 



It is nondegenerate in the sense that at no point in phase space W{y : p) 
is identically for all p. 

R2. For two states p and p' 



Part of the above requirement is that we define the appropriate measure 
dy which also fixes the constant K . The constant K = 2'Kh for continuous 
systems and K = N ioT: finite system of dimension N . This constant 
equals the volume of a "phase space cell" . With respect to this measure 
we also demand normalization condition 



Note that this is a nontrivial requirement as this implies that the left side 
of the above equation must be independent of the quantum state. 



W{y : spi + (1 - s)p2) = sW{y : pi) + (1 - s)W{y : pa) 




(12) 




1 



R3 



For any observable A on S there is a real function A{y) such that the 
expectation value (in state p) 




(13) 
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The first requirement is that W must be real. If we try to impose non- 
negativity however it becomes too stringent. As W{y : p) is convex linear on 
states in finite dimensional spaces it has a unique extension to a liner functional 
(for fixed y) on 1C{H) the set of bounded hermitian operators (observablcs). In 
infinite dimensions we need some delicate continuity arguments. Henceforth, 
we will assume linearity of W{y : p). In these specifications for distribution 
function we have not mentioned marginals. We will discuss them shortly. What 
are the characteristics of the parametric vector y? If it is to be somehow iden- 
tified with generators of classical observables its dimension must be related to 
degrees of freedom. The third item in the list gives a clue. A physical sys- 
tem, whether classical or quantum, is completely characterized by the set of 
observables O. Often O has more structure, in particular, it is an algebra. The 
difference between quantum and classical algebra of observables is that the for- 
mer is noncommutative. These algebras have minimal sets of generators. For 
example, the observable algebra of a classical system with N degrees of freedom 
is generated by generalized coordinates {qi : i = 1, . . . ,N} and the conjugate 
momenta {pi : i = 1, . . . , N}. The corresponding quantum algebra is also gener- 
ated by the operators qi and pi which do not commute. In the finite dimensional 
case we have no classical analogue. But we will be guided by this example. We 
have already discussed the close analogy between the finite-dimensional unitary 
operators X, Z and the continuous operator e*^, e"^. We show next that x and z 
are actually generators of the complex algebra of observables in the appropriate 
dimension. 

Proposition 1. Let the dimension of the system Hilbert space be N and x and 
z be as given in ([S]) and ([7]) respectively. The completion of the complex algebra 
generated by x and z equals A/„(C), the algebra of complex matrices of order N. 

Proof. The completion of algebra means that we include the limits of convergent 
sequences. In particular, X = e*^ and Z = e'^ are in the completion. We show 
that X and Z generate M„(C). Let w = e^''*/^ and Zk = uj~''Z. It is easy to 
see that 

(/ + Zfc + Z| + • • • + Z^-')/N = Dk 

where Dk{ij) = SijSjk is the diagonal matrix with 1 in the fcth row (and col- 
umn) and O's everywhere. We also see that X^ Dk = Ej+k,k where Eij are the 
elementary matrices with 1 in the ijth place and O's everywhere else. Note that 
j + k \s to be considered mod N. Thus every elementary matrix is generated 
by X and z. As the elementary matrices constitute a basis for M„(C) the proof 
is complete. □ 

We mention that the assertion of the proposition was essentially proved by 
J. Schwinger |20j in different way. The continuous quasiprobability distribution 
functions can be written as 

W{q,p)= r/(u,v)<e*^"e'P-^>e-*("i+^PMudv (14) 
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Here / is a scalar or c-number function which is usuahy interpreted as an op- 
erator ordering prescription. The Wigner distribution function is a special case 
corresponding to Weyl ordering. All this is possible because of the simple com- 
mutation properties of the observables Qi andpi. We have observed that the 
unitary operators X and Z have multiplicative relations very similar to and 
e*^ (see equation ([5])). This analogy extends further 

Z°-X^ = uj^^X^Z"; e'^^e^^^ = e'^^t^^Pe'"^ (15) 

provided a and h are integers. Of course, the second formula is valid for all real 
a and h but the first fails if both arc non-integer. This provides another reason 
to construct a discrete version of distributions functions. Henceforth we will 
restrict to finite dimensional spaces mostly. Since the operator x and z can be 
used as generators wc will assume the "phase space" spanned by y = (x, z) is 
2-dimensional. Now we can state the marginal conditions corresponding to the 
"axes". 

R4. The quasiprobability distribution function W{x^ z) have marginal distri- 
butions that coincides with probability distributions of the quantum ob- 
servables X and z: 

2 W{x, z:p) = 4,Tr(|j>0-| p) and ^ W{x, z : p) = (5.,Tr(|j> <j| p) 

X Z 

(16) 

where (resp. |j)0|) ^-^g eigenvectors of z (resp. x) with eigenvalue 

2T13/N. 



We seek a finite distribution function similar to the form p4|) above. For a state 
p in a finite quantum system of dimension N define 

JV-l 

W{x, z:p)= Yi /("^' ?i)<^"^">w~^'"'=+"^) 

m,n=0 (1'^) 

with uj = e^'^^l^ and < j, fc ^ - 1 integers 

Call the functions / in the above expression ordering functions. To compute 
the expectation values we need the following matrix elements in computational 
basis. 

(ftlA^Z-b) . (18) 
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To see the implications of the reahty condition Rl it is sufficient to verify it for 
pure states. Hence for p = |a>(a| 



N-l 



W{x,z:p)= 2 /(m,n)<a|Z-"X-"|a>a;™^ 

m , n = 



mx + nz 



m , Ti — 



2 f{N -m,N -n)i^'^''{a\X'''Z''\a)uj- 



??i,n— 1 

z : p) 



In the last step we use = = I. Since this must hold for all state vectors 
a we have 



/(7V-m,0) = /(m,0), /(0,iV-n) = /(0,n) and 
f{N -m,N- n)uj"^'' = f{m, n), l<m,n<N 

We will see later that the condition of nondegeneracy is automatically satisfied. 
Next we consider R2. Let p = Pjk \ j} {j \ k and p' = 2 P'jk b) 1 ^- Then using 

m 

Y,Wix,z:p)Wix,z:p') 

XZ 

= 2 1] /(m,n)/(m',n')pj,,+™py,,'+„,c.^»+^"'"'^-((-+'">+("+">) 

XZ jmn 

j'rn'n' ^20) 



Y^Y^f{m,n)f{N-m,N-n) 



Pj,j + mPj' .j' +N-m' 



U-j')n 



N'^ Yj Yf{'^'n)f{m,n)pj,j+^Pj',y+N-rni^ 
jmn j' 



{j-j'+m)n 



where we have used (|T9|l in the last step. According to R2 this should be equal 
to Tt{pp')/N = (Yijk PjkP'kj)/^ fo'^ choices of density matrices p and p' . This 
is possible if \ f{m, n)p is a constant independent of m and n. A straightforward 
computation yields \f{m,ri)\ = l/N"^. Setting f (171,12) = g(m,n)/N'^ we may 
write g{m,n) = a;^(™'"). We now prove existence and properties of distribution 
functions satisfying the conditions R1-R4. 

Theorem 2. For every density matrix p in an N -dimensional Hilbert space and 
UJ = e^"^^ let 

N-l 

-{mx+nz) 



W(x,z:p)= 2 /(m,n)<X™Z"V 
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Then there exist functions fim, n) such that the corresponding W satisfies R1,R2 
and R4- Moreover, for any W satisfying these conditions there are unique her- 
mitian operators a{x, z) such that following hold. 

W{x, z : p) = Tr{pa{x, z)) and p = N'^W{x, z : p)a{x, z) (21a) 

xz 

Tr{a{x, z)a{x' , z')) = —Sxx'Szz' and ^a(x, z) = I (21b) 

xz 

Given a hermitian operator T let t{x, z) = Tr(Ta(x, z)) then 

(T} = Y,Wix,z:p)tix,z) (22) 

xz 

Thus R3 is also satisfied. 

Proof. We have observed that functions f{m,n) = a;'5(™'")/iV^ satisfymg the 
relations (fT9|) provide a distribution function W{x,z : p) that satisfies conditions 
Rl and R2. To see the imphcations of condition R4 on marginals we observe 
that 

2 W{x, z:p) = Y, fim, n)<X™Z">c.-"^ ^ w"'"- 

X mn X 

= ;2 /(m, n)<X™Z" >a;-"-^„o = N^, /(O, n)<Z" ^""^ ^23) 
= NY,f{0,n)p,,(k\Z"\j}u;~-^ = Nj]f{0,n)p,,uj'-^-^^- 

jkn jn 

For the last expression to be equal to Tr(p |z){z|), the probability of finding the 
system in an cigcnstatc of z with eigenvalue 2TTZ/N, we must have /(O, n) = 
l/N"^ for all n. Computing the trace in the Fourier transformed basis |j) = f2 |j) 
we conclude that the second condition in R3 yields /(m, 0) = 1/iV^. Assuming 
^(m, n) to be a real polynomial in m and n we conclude that it must be of the 
form /3(m, n) = mna{m, n). More generally, we may take /3(m, n) = 7(to, n) + 
mna{m, n) with 7(0, n) = 7(771, 0) = 0. With this choice of (3{m, n) the first set 
of equations in (fT9|) are satisfied. The second set gives the following requirement 
on the function a. 

mn{a{m, n) - I) + {N - m){N - n)a{N - m, N - n) + 

'-f{m, n) + j{N - m, N - m) = mod TV ^^^^ 

Note that wc do not require a or 7 to be integer-valued or symmetric. There 
exist (real) functions satisfying equation p4p for all ^ m, n — 1. Simple 
solutions to these equations arc given below. 



/o(m,n)= (25) 
(J^mn — ,v2-, N even 
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where Vmn satisfies 



|fm„| = 1 and lyN-mM-ni = {-l)"'~^"i'mn (26) 



A particular choice of satisfying is 



(l-Srr^o)(l-Sr,o){m+nfN/4 



(27) 



Other choices of wih be given later when we impose more conditions on 
the distribution functions. Finally, suppose the functions / in the definition 
of W{x, z : p) satisfy the reality conditions ([T^]) and the marginal condition 
f{m,0) = /(O, n) = Then it is easy to see that for the incoherent state 

I/N 

Wi.,z:l) = l, (28) 

The distribution function is nondegenerate at each point in phase space. This 
proves the existence of solutions to equations and hence quasiprobability 
distributions satisfying R1,R2 and R4 in all finite dimensions. 

Observe that the map S(x, z) : p W{x, z : p) is real and can be uniquely 
extended to a linear map on all hermitian operators. That is, S is a lin- 
ear functional on Kh, the linear space of hermitian operators on the system 
Hilbcrt space H. Kh is a real Hilbert space with respect to the scalar prod- 
uct < A,B >= Tr{AB), A, B e Kh- Since E{x,z) is nondegenerate at each 
point there exists a unique nonzero a{x,z) e Kh such that W{x,z : p) =< 
a{x,z),p >= TT{a{x, z)p). So the first of the equations in (|21ap holds. The 
condition R2 and (|21ap together imply 

< p, p' >= Tr(pp') = 2 Wix, z : p)W^(x, z : p') 

XZ 

= W{x, z : p)Tr(a(x, z)p') = Tr((2 W{x, z : p)a{x, z))p') 

xz xz 

= < 2 W{x, z : p)a{x, z), p' > 

xz 

Since this is true for all positive definite operators p' with trace 1 we conclude 
that the second of the equation in (|21a|) must hold. Now using this expansion of 
p in the operators d(x, z) and that fact R2 again we conclude that the equations 
(j21b[) hold. Finally, to prove that condition R3 also holds observe that any 
hermitian operator T can be written in the form T = bipi — with fei, &2 > 
and pi,p2 density matrices. Then, 

<T> = Tr(pr) =< p,T>=bi< p,pi> -62 <p,P2> 

= 6iTr((2 Wix, z : p)a(x, z))pi) - &2Tr((2 Wl^x, z : p)a(x, z))p2) 

XZ XZ 

= 2 W{x, z : p)Tr{a{x, z){hpi - 62/02)) = ^ W{x, z : /9)Tr(a(x, z)T) 

XZ XZ 

= YjW{x,z : p)t{x,z) 
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□ 



The distribution functions corresponding to /q will be called (finite) Wigncr 
functions. In case of odd dimension there is one such function. But in even 
dimensions we have to be more careful in our choices. 

4.2 Explicit formulas 

The mere existence of "orthonormal" hermitian operators like a{x, z) which 
span the (real) space of observables is simply a statement about the existence of 
orthonormal bases in any Hilbert space. Two characteristics distinguish a(x, z): 
first the marginal distributions associated with them (R4) and second the way 
they were derived via the quantum Fourier transform. Our next task is to find 
explicit forms for these operators. Let 

W{x, z:p,f) = Y, f{m, n)<X™Z">w-('"-+"-) (29) 

be a quasiprobability distribution satisfying R1-R4. We have indicated explicit 
dependence on the ordering function /. From this it follows that the phase-point 
operators are given by 

a(x, z:f) = Yj f{m, n)X'"Z"a;-(™="+"^) (30) 

The fact that a(a;, z : /) form an orthonormal operator basis can be verified 
directly. The name "phase-point operator" derives from the fact that (x, z) 
may be considered as a "point" in a finite phase space. The quasiprobability 
distribution W{x, z : p, f) are simply the coefficients in the expansion of p in 
the basis {a{x, z : /)}. We will compute these operators in the "computational" 
basis = \j mod N}}, that is, the eigenbasis of the operator Z. Then X™ = 
b + s-nd = 2^ oj^" I J + m)(j|. A straightforward calculation 

then gives 

a{x, z : f)ki = (k\a{x, z : f)\l} = ^ /(^ " ^)^"^'~'^ (31) 

n 

In particular, the diagonal terms are easy. 

a{x,z: f)kk=Sk,./N (32) 

Now using the formulas p5|) for /(m, n) in the formula we get the following two 
cases for N. First for N odd, 

, ,~(k~l)x , , — (k — l)xx 

n 

Apart from ordering and normalization these are precisely the phase-point op- 
erators found in [7| for prime dimensions. Note that we do not require the di- 
mension N to be prime. If A'^ = 2r is even the calculation is a bit more involved 
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as the corresponding expression for fo{m,n) in ([25]) is not "homogeneous". We 
now have 

a(x,2 : fo)ki - Tr(|0<fc|a(x,z : /)) = c.-C^"')- j.„„c."('=+'--)/2 (34) 

n 

Evaluating these sums is not difficult but one has to be careful about the signs. 
For the choice of Vmn given in (P7|) we get 

= ) 2N^^k+i,2z k - / even 

' I ^-(fc-')^COt {T7{k + l-2z)/N)±icSc{Tv{k + l-2z)/N) _ ; qJjJ 

So we see that the quasiprobability functions given above are much more compli- 
cated in even dimension. More importantly, the phase-point operators given by 
(|33p is more sparse than the one psp for even dimension. This, in turn, implies 
that in general quasiprobability distributions are sparser in the odd dimension 
and "computationally simpler" . Let us illustrate with an example. 

Suppose a quantum circuit or protocol is supposed to produce a state |6) 
in the computational basis. Because of noise and imperfections we actually 
get a state (possibly mixed) which lies in the state space corresponding to the 
subspace K spanned by {\b ± i) : i ^ a} . From the formulas ([33|) and (|35|) it is 
easy to see that the number of nonzero entries W{x,z) in the odd case is 0{a) 
and in the even case it is O(a^). From the duality between X and Z this is 
also true if the computational basis is replaced by its Fourier transform. Since 
finding quasiprobability distribution equivalent to determining the state does it 
mean that odd dimensions are tomographically "better" ? That we should look 
at qutrits too? 



5 Heisenberg groups 

In this section we turn to our main theme: the Heisenberg groups and their 
close connections with Fourier transforms and distribution functions (see [25] 
for this connection in the continuous case) . There are families of continuous and 
discrete Heisenberg groups. Although our primary focus will be on the discrete 
Heisenberg groups we first take a look at the continuous Wigner function from a 
different perspective. We start with the (continuous) n-dimcnsional Heisenberg 
group H" whose group manifold is 1^^"+^. Using vector notation we write the 
elements as (p, q, t) where p and q are vectors in R" and t is a real number. 
The reader can easily recognize the "phase space" behind this notation. The 
group multiplication is defined by 

(p, q, t)(p', q', t') = (p + p', q + q', t + t' + (p ■ q' - q • p')/2) 

where the • denotes the usual scalar product. The symplectic structure is ap- 
parent in the above definition. By changing the parametrization of the group 
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(p, q, t) {p, q, t — pq/2) = {p' , q', t') we get the multiplication law of the (polar) 
Heisenberg group [23] . 

(Pl > 'Z'l , i i ) (P2 , 92 J ^ 2 ) = (Pl + P2 ; <?! + 92 > ^ 1 + ^ 2 + Pi 92 ) (36) 

Note that the element (0,0, i) is in the centre of the group. Let us restrict to 
n = 1 for simplicity. The Lie group is generated by the Lie algebra f)i 
with generators {p, g. A} with brackets [p, q] = A, and [A,p] = \\,q\ = 0. One 
constructs the Poisson structure on the dual space f)* in a natural way. The 
Heisenberg group plays a fundamental role in quantum mechanics. The Stone- 
von Neumann theorem asserts that the standard representation of position and 
momentum are essentially unique. In other words, the Schroedinger picture 

(p,g,t)^7(P,9,i)-e2-*e2-'(^'^^+««) 

with q4>{q') = q'i^iq') and p^{q) = ^i^-j^ is the unique representation of 
Heisenberg group under some conditions of continuity. Here ^ is the wave 
function in one dimension. Mathematically, it lives in the space H = L^(IR) of 
complex square integrablc functions (wc ignore the technical difficulties arising 
due to the unboundcdness of the operators). Since, the elements (0,0, t) are in 
the centre it is often sufficient to consider only elements of the form 7(p, q) = 
^{p,q,(S) = e'^^^''-P^^''^\ This is the reduced Heisenberg group. Let {il^a{x)} be 
a basis in H . The matrix elements in this basis are given by 

y„„'(p,g) = <^a|7(P,9)lV'o'> = <Vo|e''^^(''P+'«)|^„'> 
= Jv^4iOe2"(PP+««V„'(v)dMdi; 

These are precisely matrix elements of the Fourier transform of the phase- 
point operators in the continuous case. In particular, Vqq.(0,(7) yields Fourier 
transforms of the position probability density corresponding to the state '(/'«. 
Similarly, using the momentum representation we get the other marginal for 
Vaa'ip, 0). Since the basis was arbitrary we conclude that the Wigner function 

W{p,q) = J<7(u,w)>e-2"(P"+9")dudw 

when integrated over the strip between q = ci and q = C2 gives the probability of 
the particle in (pure) state "0 to have its position between ci and C2. Explicitly, 

f 'dq W{p,q:m>lp 
Jci J — x 

yields the probability that the position observable has value between ci and C2 
and similarly for the momentum observable. This is easily seen by expanding 
1^) in the position basis. What do we get if we integrate over an arbitrary 
strip, not necessarily parallel to the p or q axes, say the lines ap + bq = ci and 
ap + 6g = C2? The answer is well-known and is discussed in [7] and [32j . But 
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we look at it from a different perspective. First, put ap + bq = p'. This defines 
a family of parallel lines p' = c in the p-q plane. Another line cp + dq = q' 
does not belong to this family if and only if ad — be ^ 0. Thus, the matrix 

C — invertible and defines a change of coordinate in the phase plane. 

Then the form up + vq = u'p' + v'q' where (u', v') = (u, v)(^. This in turn defines 
a transformation on the Lie algebra generated by p, q: 

7(u,w) = e2"("P+"«) = e2^*("'p'+»"3')^ " ^ 

If the transformation p ^ p' ,q ^ q' were an automorphism then the p' and q' 
have the same commutation relation as p and q. This will happen if and only if 
det C = ad— 5c = 1. But then if we change the variable of integration to p', q' the 
measure remains unchanged (| det C| = 1). We can now carry over the argument 
from the case of axes marginals and conclude that integration of W{p, q : 
over a strip between ap + bq = c\ and ap + bq = C2 gives the probability that 
the observable p' = ap + bq will have value lying between ci and C2. Let us 
observe that ( e SL{2,M.) = Sp{l,M.) where SL{n,M.) is the group oi n x n real 
matrices with determinant 1 and Spin, M) is the real symplectic group of order 
n. In general for Sp{n, M) is a subgroup of the automorphism group of H" and 
is different from SL{2n,M.). 

Now we turn to the discrete Hcisenberg group H. We define a presentation 
of the group in terms of generators and defining relations [5]. H is generated 
by {x, z,7}. The defining relations are 

zx = 7XZ, 7X = X7 and 7Z = Z7 (37) 

The advantage of this approach is that any map from the generators of a group 
H to another group K which satisfies the same defining relations as above can 
be uniquely extended to a group homomorphisms H ^ A'. A simple realization 
of the group over integers is given by the set 7? . The multiplication is defined 

by 

(ji,fci,ii)(j2, ^2,^2) = O'l + J2,fci + k2,ti +t2+ iifca) 

The generators are z = (1, 0, 0), x = (0, 1, 0), and 7 = (0, 0, 1). If we specialize 
to Zjv, the integers modulo N we get corresponding Hcisenberg group Hjy with 
generators X, Z, and 7 and the relations 

= = -f^ = e (identity), 7X = X7, 7Z = Z7 and ZX = -/XZ (38) 

Since H^r is a finite group its finite-dimensional representations are unitary 
and completely reducible. Let be a representation of H or Hjv on a vector 
space V of finite dimension. Wc say that the central element 7 acts maximally if 
the order of 0(7) is dim(T^). The following theorem characterizes representation 
of H (H^v) and their relation to QFT. 

Theorem 3. Let (p be a (unitary) irreducible representation o/H^r on a finite- 
dimensional space V . Let r be the automorphism ofUN (and H ) given by X ^ 



21 



Z, Z ^ X and ^ — > 7 ^ and (j>' the representation defined by (l>'{g) = (ftirg). 
Then the following statements are true. 

1. 4){'~i) = u), a primitive Nth root of 1 and the eigenvalues of (t>{Z) and 
(j>{X) are {w*'' : ^ iV — 1}. 7 acts maximally if and only if (/) is faithful 
(one-to-one). 

2. <j> and (/)' are unitarily equivalent: (f)' = ilipft^ and ft is the quantum Fourier 
operator. 

3. Any unitary irreducible representation ip of the full discrete Heisenberg 
group H in V in which the order of 7, 0(7) = dim (V) = K is equivalent 
to an irreducible faithful representation ofUx. 

Proof. Assume first that 7 acts maximally. Since (j) is irreducible and 7 is in 
the centre it must act as a constant (Schur's lemma). As order of 7 is A^, 
7 = Lol where w is a primitive A'^th root of unity. Since H^v is finite we may 
assume the representations to be unitary. Let a be an eigenvector of (j){Z) 
with eigenvalue c. As = e, c must be an A^th root of 1. Consider the set 
S = {a, (p{X)a, (t>{X^~^)a}. As 

= (l){-f'')(j>{X'')<j>{)a = cuj''(l>{X'')a 

(j){X'^)a = (j){X)'''a, k = 0,1, . . . , N — 1 are eigenvector of 4){Z) with eigenvalue 
cw'^ . These eigenvalues are distinct roots of 1 and hence S is linearly independent 
and a basis of V. We can reason similarly for 4>{X). The converse is trivial. If 
7*^ = / for k < N then (p cannot be faithful. 

Next we recall some facts from the theory of characters associated with 
representation of a group [33| . If p is a representations of a finite group G 
on a finite-dimensional vector space V, the character Xp is a scalar function 
on G defined by Xpid) = Tr(p((7)). It is constant on conjugacy classes. If we 
have two characters Xp ^nd Xp' corresponding to representations p and p' then 
their scalar product is defined as {Xp,Xp') = (V^) Sgeo Xp(5)Xp'(5)- It is a 
fundamental result that two irreducible representations p and p' are (unitarily) 
equivalent if and only if (xp, Xp') ^ 0. We apply this to the representation </> and 
(j)' of Hat. First, observe that since ZX'^Z^Z-^ = 7™X'^Z", X0(^'"^") = 
<j-'^X4>i^^ Z^) which is possible iff either m = or X0(X"'Z") = 0. Conjugating 
with X we conclude that Xtj> is nonzero only on the centre of Hjy. Hence, to 
prove equivalence of (f> and (/)' it suffices to show the scalar product of X4> and 
X4>' is non-zero. But (j) and 4>' have the same effect on the center (generated by 

7) of Hat. Hence (Xp,Xp') = (1/^) Sfc X0(7'')X0'(7'') = 1- Since (j) and (j)' are 
equivalent there exists a unitary map : V ^ V such that <f>'{g) = il^<f>{g)il. In 
particular, (/)'(Z) = 0(rZ) = ^(X) = r2V(^)^^- Now let {|j> : j = 0, . . . , } be a 
complete set of eigenvectors of (t>{Z) with (j){Z) |j> = uj^ |j> and similarly let {|j)} 
be an eigenbasis of (j}{X). Then 0(Z) = |j><j| and (j){X) = J]^u;^ 

Observing that : j, fc = 0, . . . , A^ — 1} form a basis the space of operators 

on V it is easy to check that = | j) {j \ . We have also seen that <j>{X) | j) = 
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\j + 1). From these and the normahzation (j|0) = 1 we get (j|ri|fc) = j ^pn. 
We have proved item 2. 

To prove the last assertion we again start with an eigenvector a with eigen- 
value a of Note that we can no longer assume that a is an KXh root of 1. 
However, the hypotheses that order of ^/'(7) is K implies that a, V'(x)q:, . . . , ?/'(x^^~^)a 
are eigenvectors of V'(z) with distinct eigenvalues aip{'yy , j = 1, . . . , K—1. They 
must be then independent. This implies x^a = a and hence = 1. Hence the 
eigenvalues of x must be Kih roots of 1. Interchanging the role of x and z we 
conclude that a must be a primitive Kth root of 1 and the assertion follows. □ 

Note that the condition on 7 (maximal order) is necessary in case of the group 
H and Hat. For example, let iV = 8,^(7) = -/, p{Z) |0> = |0>,p(Z) |1> = - |1>, 
p{Z) 1 2) = 1 2) and p{X) the cyclic permutation. Then p is an irreducible repre- 
sentation of Hg . We see the connection between representations of the Hcisen- 
berg groups and the QFT. For a vector a wc write a = U,a for its Fourier trans- 
form. The Plancherel formula jjap = lap is simply stating that the Fourier 
operator Q, is unitary. We also note that since the representations of the group 
H is equivalent to Hat when 7 acts maximally it will be sufficient to consider 
H^v in a fixed representation space. However, when we are dealing with differ- 
ent representations (for example, taking tensor products) we have to deal with 
the full Heisenberg group. The condition that 0(7) = dim(y) is special case of 
general irreducible representations of H. It is sufficient for our purposes and 
will be implicitly assumed. Henceforth, for a fixed representation p we simply 
write the action of a group clement g as ga instead of p{g)a if the context is 
clear. 

Now wc turn our attention to distribution functions. We have seen that the 
distribution functions can be given an alternative characterization in terms of 
phase-point operators. The formula ([30|) for these operators implies that they 
are linear combination of the unitary operators of the group. Thus, we look for 
them in the group algebra. Recall that for a group G the group algebra C{G) over 
complex numbers is the set of formal finite linear combinations Ci^i, gi e G 
and Cj 6 C. The algebra product is defined as 

i j ij 

Any representation of the group is a representation of the group algebra and vice 
versa. Now for a unitary representation p of G on a finite-dimensional vector 
space the character Xp of the representation induces a scalar product on C (G) . 
Thus 

(m> i^) = Xp(M*i') where = ^ Qg^, = ^ c'^g'j and fi* = ^ Qg,^^ (39) 

This is indeed a scalar product on C{G). The resulting norm coincides with the 
Hilbcrt-Schmidt norm on the corresponding operators on V . Call an element 
/i 6 G(G) self-adjoint if /i* = ji. Let G = Hat or H. Since the central element 
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7 acts as a scalar we write elements of C{G) in the form „ Cm„X'"Z". For 
a representation (jj of H^r with 7 acting as ujI consider the following elements 

A{x, z) = Y, c™„w-'"^+"^X™Z", X, z = 0, . . . , iV - 1 in C(Hjv) 

mn 

We demand that the set Q = {A{x, z)} be mutually orthogonal, self-adjoint and 
satisfy the following: the elements P{x) = ^(^' ^) Qi^) = l^x ^{^: ^) 
are projections, that is, P{xY = P{x) and P{z)^ = P{z). We call such a set of 
elements in C(G) a Wigner set. We have the following theorem. 

Theorem 4. For a representation (j) o/H^v on N -dimensional space V , Wigner 
sets exist in C(Hjv)- If G is a Wigner set then for A{x,z) e G, A{x,z)/N are 
phase-point operators. In other words, given a quantum state p the function 
W{x,z : p) = li:{(j3{A{x, z))p)/ N is a distribution function. Conversely given a 
distribution function W{x, z : p) on V there is a unique Wigner set A{x, z) in 
C(Hjv) such that W{x, z : p) = 1t{4>{A{x, z))p)/N . Wigner sets are translation 
invariant in the sense that the transformation Cmn ~> Cmnw"'"^^", a,b e Z]y 
permutes the operators A(a;, z) in a Wigner set. 

Proof. The proof is similar to that of Theorem [2j We only sketch some of the 
basic arguments since we are dealing with group algebras. First, the self-adjoint 
property implies conditions like (|19p with Cmn in place of /(m, n) since the 
are independent in C{IIn). Let us compute the scalar product of two 
elements from G- Assuming now self-adjointness we have 

iAix',z'),Aix,z)) = 



Tr 



\ ^ — mn "v^m-fm y-n + n , — {fn x +n z ) —{mx + nz) 

/ J ^m'n'^mn^ Zj Lu UJ 

m,n 

\m' ,n' / 



m.n 
\m' .n' 



rp I , — m{x—x')—n(z — z') 

\2^-m{x-x')-n(z-z') 

In deriving the second step we use the fact that Tr(0(X^ Z'^)) = unless j = k = 
mod N. The last expression will be proportional to Sxx'Szz' if |cmnP = a 
constant. Wc will fix K shortly. Hence, we assume that c,nn = Kuj^"^'\ For the 
last requirement we have 



^ z(m+n) 



This would be possible if con = K = for all n. We have already proved the 
existence of functions satisfying these conditions in Theorem [5J The fact that 
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W{x, z : p) = Tr(A(x, z)p) is real follows from sclf-adjointness. The orthogonal- 
ity property implies R2 in Section 14.11 

Tripp') =Y,W{x,z : p)W{x,z : p'). 

xz 

Finally, the property about marginals is equivalent to showing that P{z) and 
P{x) represent projections on |z> and |i> respectively. We can deduce this 
directly from the fact that 

0(Z) =Y,^^ and </>(X) =Y,^^ 

3 3 

The proof of the converse is straightforward. 

To prove the last statement let c'„„ = c„i„w°™+''". Then 

mn 

\ ' „ , —in(x—a]-Yn{z — h) ryn 
mn 

= A(x' , z') where x' = x — a and z' = z — b 

The assertion follows from this and the proof is complete. □ 

We see the correspondence between orthogonal sets in the group algebra 
C(H7v) and distribution functions. We have seen that the QFT arises out of a 
particular automorphism t of the Heisenberg group. So we expect the general 
automorphisms of H^v and H to contain more structure and information relating 
to QFT and distribution functions. It is easy to see that any two representations 
of Hjv in which 7 acts maximally and has the same value are equivalent. In 
particular, if a is an automorphism of H^r that fixes 7 and cj) is an arbitrary 
representation then cf) and • a arc equivalent. If {A(x, z) : ^ x, z ^ N — 1} 
is a Wigner set then {aA{x,z)} is also a Wigner set. So we can generate new 
Wigner sets by automorphisms. As the value of cr on X and Z determines it on 
C(HAr) let a{X) = X^Z^ and a{Z) = X'^Z'^. We must have 

(^a^6)7V ^ ^a6iV(7V-l)/2 ^ (^c^d)JV ^ ^cdA'(JV-l)/2 ^ ^ 

^c^d^a^b = 7'^'^-f«=X''Z^X=Z'^ = 7X"Z^X=Z'^ 
for CT to be an automorphism. The second condition implies that the matrix 

«'•[: % 

has determinant 1. That is, e SL{2, Z^), the set of matrices with entries 
in the ring Z^r and determinant 1 mod N. Conversely, given M e SL(2, Zjsr), 
N odd we can define an automorphism cm as above. If N is even this simple 
definition of aM does not work in general. For example, if ab is odd then 
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{X°'Z^)^ = 7^/2 = —1. This is reminiscent of half-integral representation of 
rotation group {SU{2) actually). Hence for even N we have an automorphism of 
H rather than H^v- Note that in this case for any M e SL{2, Z^), {aM{g))'^^ = 
1, g e Hat. There is however a proper subgroup of SL(2, Zjq) which induces an 
automorphism of Hat. Alternatively, for even N define the function 



sgn(u) 



0, w e Z and u even 

1, u odd 



Now we can define the automorphism on H^v, Me SL{2, Z^) for odd N 
and for even N we define it on the representation space. 

N odd 

^^^^^^ - |^sg„(a.)/2^a^.^ ^ (40) 

Let A{x^ z) be a Wigner set and a representation of H^r. Then we have seen 
that a distribution function is defined by 

W{x,z : p) = Tr(0(A(a;,z))p)/iV 

This means that if and 0' are equivalent representations connected by a a 
unitary operator U and W and W are the corresponding distribution functions 
then 

W{x, z: p) = W\x, z : pU) (41) 

In particular, if M e SL{2, Zf^f) then it induces an equivalent representation (/jm- 
In case A'^ is odd (pM is the representation that is given via the automorphism 
generated by M. In even dimension <j>M is defined by (|40)) above. Now let us 
look at other "marginals" of a distribution function W{x,z : p). One way of 
constructing such marginals is via a finite Radon transform [571 HH] ■ Thus for 
/ : Zjv X Zat ^ C define "lines" 

Sab = {{x, z) e Zat X Zjv : ax + hz = Q mod N, gcd (a, h,N) = l] 
W{x\z':p)= W{x,z,:p) (42) 

The condition gcd (a, &, N) = 1 ensures that the "line" ax + bz = t has a solution 
for all t 6 Zm- The "coordinate axes", for example, correspond to the sets S'lo 
and S'oi. The function W{z') in is a Radon transform of the distribution 
function Vl^(a;, z) and each pair (a, 6) e Zat x Zat such that gcd (a, 6) is invertible 
in Zjv, defines such a transform. Let x' = ax + 6z. Let (c, d) e Zat x Zat such 
that the matrix 

a h 



M=\^^ e5L(2,ZAr) 

and set z' = cx + dz. In the following it will be convenient to use vector 
notation. Thus ^ = (x, z) e Zjv x Zat is a two-dimensional "vector" 0. We will 



^To be accurate Zat x Zat is a module over the ring Zf^. 
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also occasionally use the component notation: ^ = (^1,^2)^- So the distribution 
functions and phase-point operators may be written as : p) and A{() 

respectively. Then a marginal with respect to the second component is given by 

W{z' ■.p)= ^(^^^ ■ p) 

Here z' = M^2- In analogy with the continuous case we require that W{z' : p) 
be a probability distribution with respect to z' . More precisely, in the repre- 
sentation ipM-^ of the Heisenberg group corresponding to the automorphism 
induced by M'^, W{z'; p) gives probability distribution of the quantum observ- 
able — iln0j\f-i(Z) in the state p. However, there is a sharp difference between 
the distribution functions in even and odd dimensions. The general marginal 
condition holds in odd dimensions for the Wigncr distribution function defined 
by (j25p but not for even dimensions. For even dimension we have more com- 
plicated formulas for the marginals of the Wigncr function. In fact wc will 
show that in this case no distribution function satisfying conditions R1-R4 in 
subsection 14 . 1 1 will satisfy the general marginal condition for all AI e SL{2, Z^). 

Theorem 5. Let (p he a representation ofH^ on V. Let 

Aix, ^) = ^ E /("^, (43) 

m,n 

be a Wigner set and W{x, z : p) = Tr{A{x, z)p) be the corresponding distribution 
function. 

M=(^l ^^eSL{2,Z^) 

We say that simple marginal condition ( with respect to M ) is satisfied if the 
marginals 

W{z':p)^ 2 W{M-'^:p) (44) 

are probability distribution in the eigcnbasis of the operator (pM^Z) = (f)(aMiZ)) = 
(j){X'^Z'^). Then the following statements hold. 

1. The (simple) marginal condition is satisfied if and only if 

Am{x,z)^A{M-^Ci 
is a Wigner set for the representation (f>i\j. 

2. If marginal conditions are satisfied for all M e SL(2, Zfi) then the Wigner 
sets (and distribution functions) are determined uniquely up to transla- 
tions. 

3. If the dimension N > 2 is even it is not possible to satisfy the marginal 
condition for all AI e SL(2, Zn). 
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4- In odd dimension Wigner functions are given by 

W{X,Z) = 2 ^™"(JV+l)/2^m^n^-(,nx + nz) ^45) 
m,n 

up to translations. In that case W{z' : p) = {az' -vcdl2\p\'^z' -^cdj"!) where 
\oLj) are the eigenvectors o/crj\/(Z). 

Proof. Since gcd (a, b, N) = 1 there exist integers c, d and k such that ad—bc + 
kN = 1 and hence ad - be = 1 mod N. For a matrix M let M' = (M^)-^ 
The first assertion in the hst is relatively straightforward. If the equation (|44p 
is satisfied then Wm{x,z) = W{M~^x, z) is a distribution function. The 
condition Rl (reality) is clear, R2 follows from the state transformation equa- 
tion (|4T|) and (|44p gives the marginal condition R4. From the correspondence 
between distribution functions and Wigner sets the first assertion is clear. 

Suppose the marginal conditions are satisfied for some M e SL{2, Zn) given 
above. Using the formula 

^j^ti^i'-jrn ^uvra{rn—l) 12 j^ura ^vm (46) 

we get setting ^ = (m, n)^ e x Zm 

= TiTi f{m,n)X"'Z"uj-^^'^-^ = E E /(Af^C')^™^"^"^''^ 

m,n MCi C'=M'C 

aM(X)'VM(Z)"^-'™"+""\ (a = x) 
= 2 /(cn, dn)w-(^'^"("-^)/2)^^^(^)n^-«. 

n 

= 2g(n)Z'", 5(n) = /(cn, dn)^-^"''"^"-!)/^) ^^j^^j ^ aM{Z)uJ-' 

n 

We require that the operator T = "^^^ gin)Z'" be a projection: it must be 
hermitian and satisfy = T. Hence, we must have 

2 g(m)g(n)Z'"+™ = "^)^" = 

Since Z' like Z has no repeated eigenvalue its minimum polynomial is the char- 
acteristic polynomial — 1. Hence the operators I, Z' , Z'^ , . . . , Z'^^^ are 
linearly independent and we must have 

Yi9{'m-)9{k - m) = g{k) 

m 
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But the left side is the convolution of the function g with itself. Taking (finite) 
Fourier transform of both sides we get = g. Then g{m) = 1 or 0. This 
implies that 

i 

where ti e Zj^ are the values at which g = 1. But from condition R2 we infer 
that \g{x)\ = 1 and since Rl implies that g{0) = 1 we conclude that there must 
be exactly one term in the above sum: 

Putting n = 1 this implies that f{c,d) = uj*^''''''-^ whenever gcd(c, d, A^) = 1. 
Hence we rewrite the above equation as 

/(cn, dn) = f{c, d)"t^c'i«("-i)/2 and so f{cn, n) = f{c, l)"t^="("-i)/2 (43) 

Now suppose N is even. By the definition of Wigncr sets they must be indepen- 
dent since the operators are mutually orthogonal. Consequently, the function 
/ must be periodic with period N and since /(c, 1) = a;*^^'^\ t{c,l) must be 
an integer. Putting n = iV in the second equation in (pS)) and noting that 
f{x,0) = 1, Vx 6 Z we get a contradiction when c is odd for the right side is 
^1. Hence it is not possible to have Wigner sets satisfying all simple marginal 
conditions. 

Next suppose that N is odd. Then 2 has an inverse (iV + 1)/2 in Zn- It is an 
easy verification that the function /(m, n) = cj™"(^+^)/^ satisfies the functional 
relation (|19p . To prove uniqueness we assume that t(m,n) can be extended to 
all Z and that it can be expressed as a polynomial in m and n with integer 
coefficients (which may depend upon N). Since /(m,0) = /(n,0) = 1 we may 
assume that the polynomial is of the form t{m, n) = ^;™"[''o+s(m,n)] where aq is 
a constant and g{'m, n) is a polynomial without constant term. Then we have 

f{cn,n) = ^™^(oo+3(cri,Ti)) ^ ^ric(ao+g(c,l))^cri(n-l)/2 

Since this must be satisfied for all n we must have oq = (N + l)/2 and g = 
mod N. This proves uniqueness up to linear terms. 

The last statement is easily derived from the above proof of existence and 
uniqueness of distribution function satisfying all the marginal conditions for odd 
N. □ 

We note that a similar relation holds for the marginal distribution over x 
when wc average over the variable z. In fact, satisfaction of marginal conditions 
under the full SL(2, Zn) for one variable implies the same for other. In even 
dimensions there exist no distribution function satisfying all marginal condi- 
tions. Therefore, we have to relax some of the conditions of the theorem to 
get the marginal distributions. Let us recall why the marginal conditions are 
desirable. One of the main reasons is that by determining sufficient number 
of marginal distributions we can reconstruct the state. The simple marginal 
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condition stated in the Theorem [S] is satisfied (see (|33]) and the statement that 
follows it) the marginal distribution corresponds to probabilities for a complete 
projective measurement in a suitable basis. In even dimension we have three 
options. 

1. Wc do not require that the Wigner set be independent. Then the repre- 
sentation of the Hcisenbcrg group Hat need not be irreducible. This was 
the approach adopted in }25| . 

2. We drop the conditions that the marginals be simple type. As will be 
shown next we can still determine the "marginal" distributions from the 
measurement probabilities. 

3. We do not demand that the marginal condition be satisfied for the full 
SL{2, Zn) but only for a subset. We show that in case N = 2^ there is 
such a subset and the marginal distribution for it are sufficient to recon- 
struct the distribution function. 

Wc start with the first option [25l [16]. Since the operators A{x,z) are no 
longer independent the function / (as a function on Z) is not required to be 
periodic and the labels (x, z) can take any integer values. A minimal extension 
is obtained by looking at the basic recurrence relations pS)) . The problematic 
factor a;'="("~i)/2 is periodic with a period 2A'' (as function on Z). The same 
relations then suggest that we take f{m,n) = w"'"/2 where w^/^ is a primitive 
27Vth root of 1. Hermiticity of phase-point operators then require that we now 
define them as 



Because of redundancy these operators are not uniquely determined (up to linear 
factors). But we can modify the proof in Theorem [5] for odd dimension to 
determine the possible solutions in this case. 

Next we look at option 2. We defined a family of distribution functions say 
W{x, z : p,v) in the even case in (j25p depending on some function v. The func- 
tion v is arbitrary apart from the condition ([25]) . Let Wq denote the special case 
when V is given by (p7|) . Of course, Wq does not satisfy the marginal condition 
but the results below show how it may be computed from the measurement 
probabilities. 

Proposition 2. Let V be an irreducible representation space ofJi]\[, N even. 
Let 



Let u = gcd(t, N) where t = c if c even and d otherwise. Suppose N/u is even. 



A{x,z)= 2 



OJ' 



nin/2 \rm ryn 
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Then 



W{z':p,v)= W{M-^x' ,M-^z' : p) 

x'eZ„ 



(49) 



where \x\ is the greatest integer ^ x. In particular, for Wq with v^n given by 
(123) we have 

Woiz' : p) = 

{(a^, cd-sgn(cd) \p\a^, cd-sgn(cd) ) + jy , cd-5gn(cd) \ p\Oi N_ , ^, _ cd-ssn{cd)_}) /2 

" 2 2 2 ^ 2 2 2 

j n odd 



where h{i, z' 




N 



if cd odd 



. 27r(£i^+j-z') .J, , 

sm — ^ — ^ if cd even 



(50) 



The proof is given in the appendix. Observe that the if TV = 2^^, K > 1 
then A^i = N/u is always even. We can also write the appropriate formulas 
for the case A^i odd. We avoid doing so as they are even more complicated. 
We can also simplify the trigonometric sums in (|49p . However, note that if 
we know the probabilities then in principle the Radon transform W{z : 

p, M) can be computed by evaluating these sums. In case of odd dimension 
the expressions for the marginals arc simpler but we still have to estimate the 
probability distribution in the basis {jaj)} defined above. If we have these the 
probabilities doing the sums in the even case is routine. So, is there a deeper 
reason for imposing the marginal conditions on the distribution function? Two 
possible reasons could be simplicity and some theoretical insight. 

We consider the second option listed above for dimension N = 2^ only. 
Thus we aim to construct a distribution function which satisfies the marginal 
conditions for only a subset of SL{2, Z^). The theorem below gives an explicit 
formula for this important case. Thus let Li c SL{2, Zn) be a subset consisting 
of the following matrices. If M e Li then each row has at least one entry = 1 
and if the diagonal entry is 7^ 1 it is even. 
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Theorem 6. Let TV = 2'=. Define 



even 



m>0 
n odd 



2 (-1) I. '""w l^j^m ^ny^mn/2-{mx + nz)'^ 



n even >0 



(51) 



where the expressions like (rnn~^) are first computed modulo N in the residue 
class {0, . . . , A'^— 1} and then treated as an integer. [x\ denotes the largest integer 
less than or equal to x. Then Wi satisfies the conditions R1-R4 and for every 
M e Li, Wi satisfies a simple marginal condition with respect to the variable x: 

Wi{z':p,v)^ Wi{M-^x',M-^z' : p) 

_ |(«./ c + .gn(c) c + .gn(c) ), d=l ^"^^^ 

\(ct^j_d \p\a^,_dy , c = 1 and d even 

Proof. Wc first note that the notation n^^ makes sense in the ring Zf^ since 
every odd n is invertible. The reality condition Ri is seen from the following 
simple observation. For Q ^ m,n < N let mnr^ = kiN + ni, ki ^ and 
< ni < N. Since {N — m){N — n)~^ = {mn~^) mod N we get 

{{N - m){N - ny^){N - n) = (to - 1 - ki)N + m 

This implies that ki has same (opposite) parity as y{{N — m){N — n)^^){N — n)\ 
if TO is odd (even). Hence 

We can argue similarly for to odd and n even. Hence the reality condition (|19p 
is satisfied. It is clear that Wi is normalized. The other conditions easily follow 
from the definition and the analysis of these conditions in Section 14.11 Finally, 
the simple marginal condition with respect to x is seen to be satisfied as follows. 
From Theorem [5] and Proposition [2] we note that we have to consider pairs of 
the form (cti, dn), where (c, d) is the second row of Af , in the calculation of the 
marginals. Using the notation of Proposition [2] we set 



1 TO, n even 

I (mn^^)n I j j 

Tv 



(nm I j i 

^ ' — J TO odd , n even 
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As the matrices belong to Li we consider two cases. If the diagonal element 
d = 1 then the only terms in the sum yielding Wi that contribute to the marginal 
are indexed by ((cn), n) where n runs through Zj^ and (cn) is calculated mod N . 
The case c = is already covered. If c 7^ then from 



Wi{z':p,iy)= Wi{M-'^x' ,M-^z' : p) 
j " 

j n odd rx even 

= 2p^.^.^^((c+sgn(c))/2+j-.')n = _ C + Sgn(c) _ , _ C + Sgn(c)^ 
j n 2 2 

For the case c = 1 and d even the terms in which n is odd drop out from the 
sum for Wi{z' : p, v) and the proof is similar to the first case. □ 

We note that the subset Li of matrices from SL{2, Zjsi) cannot be extended 
arbitrarily preserving the property of simple marginals. For example, if we 
admit matrices with c = 1 and d odd then we get a factor of sgn([ *-'^jy'*" J) 
instead of sgn([^J). The two need not be equal. However, as we will see below 
the set Li is sufficient to determine Wi. 



5.1 Inverse Radon transform and state determination 

In the previous section we have seen that the finite Wigner distribution function 
enjoys a rich variety of marginal properties. We can use this to determine the 
former. This is equivalent to inverting a finite set of Radon transforms. From 
the distribution function we can determine the state. The invertibility of the 
Radon transforms also shows that the Wigner distribution function is unique 
up to a translation. In the rest of the section W{x, z) will denote the Wigner 
distribution function. Replacing the matrix by 




dct M = 1 mod N 



in Theorem [5] we rewrite the basic Radon property stated in (|44|) and the state- 
ment following. For example, for odd N 

W{z : p, M) = 2 W{Mx, Mz : p) = Tr(|a,_ac><«3-ac| p) (53) 

The problem is to reconstruct W{x^z) from W{z : p,M). Call the later 
the Radon transform of W with respect to the matrix M. The idea is that 
W{z : p,M) is the probability distribution of the observable —i\n{X~'^Z°') in 
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the odd case. In case of even dimensions it can be computed from the distri- 
butions. Assuming that these distributions can be approximately determined 
experimentally we can reconstruct W and hence p. We have seen that in odd 
dimension N there is a distribution function satisfying simple marginal condi- 
tions for every M e SL{2, Zjsi) and in dimension N = 2^ we have only a subset 
of SL(2, Zff) with simple marginal conditions. We give explicit formulas for 
these two cases. First some notation. For a subset S of some set let xs denote 
the indicator function: xs(^) = 1 if a; e S" and otherwise. In the rest of the 
section we use the boldface vector notation to denote a member of Zjq x Zm 
and other non-bold letters to denote "scalars" belonging to Zm- For example, 

^l= [^^] , /ii,/i2 6 Zat 



Given M e SL{2, Z^) let Ci(M), i = 1, 2 denote the column vectors of M . Let 

S,{M) = {C,{M)x : X e Zn} c Zat x Zjv, i = 1, 2 

Theorem 7. Any distribution function W{x, z) can be uniquely determined 
from the (finite) set of Radon transforms W{z : p,M) where 

M=(^l ^^eSL{2,ZM) 

In particular, for odd dimensions the Wigner function given in (j45p and any M 
we have 

W{ct,-at) = ^2<a,_„,/2(M)|p|«,_,,/2(M)>c^^* (54) 

Z 

and for N = 2^ and M <^ Lx 

(^<^az'_|.(Af)|p|a2'_|(Af)^cj^*, c= —1 and a even 

where W is the Fourier transform of W in Zj^ x Zjv and |aj(M)) are the 
eigenvectors of aM-^{Z) = X~'^Z'^^ In either of the cases the Wigner function 
W or Wi can be reconstructed from the marginal distributions. 

Proof. Write the Radon transform ((53|) as 

Wiz ■.p,M)= Yj W{Mx, Mz : p)= ^ W^(C2(M)z - x : p) 

where 5*^ = — 5*1 We can write this as a convolution. Thus 
W{u : p,M) = xs.(u)2w^(u-x : p)Xs[(ai){^) 

X 

= Xs.(u)(VF*Xs;(A/)(u))> = {C2{M)z} 

^Recall that here M replaces M^^ of theorem [5] 
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Now we take the finite Fourier transforra of the above equation in the group 
Zf^ X Zm [22]. Recall that the Fourier transform of a complex function /(u) on 
Zjv X Zn by / is a function on the dual group {Z^ x Zat)*: 

Ml ,M2 

Using the fact that W = Wxs^ and that the Fourier transform of a convolution 
is a product and vice versa we have (suppressing p and M) 

Ml ,U2 

{u-.aui -{-01-2—0] 

t 

where F(t) = W{ct^—at) and F is its Fourier transform in Zn- In proving 
the above we use the following facts: Xs;(J\/)(i^) ^ iff avi + cv2 = and the 
solution to the congruence equation avi + cv2 = mod N is given by the set 
{(ct,— at) : t G Zj^. This follows from a similar result for linear Diophantine 
equations [31] and the fact that gcd(a, c,N) = 1. We also use det M = ad— be = 
1 in the las but one step. The factor ^/N appears because of the normalization 
used in our definition of FFT. It now follows that 

Wict,-at) = ^Y^Wiz),,^' (56) 

z 

This formula is valid for any distribution function. Let now N be odd. Com- 
bining this with the equation (j44p in Theorem [5] we get (j54p . Similarly, when 
TV = 2^^ and e Li wc obtain ([55]) . Note that the formulas in ([5^ arc valid 
under the assumption that M e Li (see the footnote 2] above) . 

Wc next show that it is always possible to find a,c e Zn such that gcd (a, c, N) = 
1 and the "lines" {{ct, —at) : t e Zj^} cover the "plane" Zjsr x Zjy in the two cases 
above. When TV is odd this is obvious. If iV = 2'^ consider {x, y) e Z^ x Z^. For 
< j < A'^ let hj denote the highest power of 2 that divides j, that is, j'/2''j is an 
odd integer. If ^ hy then we put a = 1 and c = — 2''=="''" (y/2''»)~-^ where the 
inverse is evaluated in Zjy and we assume that y ¥= 0. Then {x,y) = (ct, —at) 
for t = —y. If hrc < hy then put c = 1 and a = —2'^y~^''{x/hx)~^- We have 
therefore shown that in alljiiese cases the Radon transforms together can be in- 
verted for from the values W{ni, ^2) so obtained we can take the inverse Fourier 
transform and the last assertion of the theorem is proved. □ 
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We can thus recover any distribution function W{x, z : p) and consequently 
the state p from the Radon transform data which are in turn probabihty distri- 
bution of measurement in appropriate bases (see The theorem shows the 
existence of an inverse transform corresponding to the set of Radon transforms 
of W, each corresponding to an element M in the group SL{2, Z^). But we do 
not need all the Radon transforms. What is an optimal subset Q cz SL{2, Zm) 
that suffices to determine the state uniquely from probability distributions cor- 
responding to measurements in appropriate bases? This question can only be 
satisfactorily answered in the context of prior information about the state. One 
can show that without any such information the cardinality of Q is 0{N). Even 
then we have a lot of freedom. We can use our choices so as to ensure optimal 
measurement. Recall from Theorem [5] that the Radon transforms are given by 
probability distribution (corresponding to a state p) in the basis that diago- 
nalizcs the unitary operator X'^Z'^. The only condition imposed on the pair 
(c, d) e Zjv X Zjv is that gcA {c,d,N) = 1. We can often compute this basis 
explicitly. Then we can use quantum circuits to transform our original "compu- 
tational basis" to the required basis. A criterion for the choice of (c, d) could be 
those that minimizes the size of the circuit. For example, if = 6, the choice 
c = 3, d = 2 leads to a particularly simple basis. The analysis becomes simpler 
if the dimension iV is a prime power. We aim to address these issues in future. 

5.2 Distribution functions and quantum information 

In this section wc discuss some potential applications of distribution functions 
in quantum information processing (QIP). This is a developing area and wc only 
sketch how our formalism may prove useful in various areas in QIP. For this it 
is best to view the distribution function as coefficients in the expansion of the 
state in some orthonormal basis in the space of operators, in particular, the basis 
consisting of phase-point operators. First we generalize to automorphism groups 
of the group algebra C{H.m)'- a linear isomorphism T : C(Hjv) — > C(Hjv) such 
that T{xy) = T{x)T{y) is bijective. It is sufficient to check the last condition 
for the generators X, Z and 7. We will consider only those automorphisms for 
which T{j) = 7. Then T{X), T{Z) and 7 generate a group isomorphic to H^v 
provided T{X)^ = T{Z)^ = 1. In particular if c e C{IIn) is invcrtible then 
the map T{x) = cxc~^ is an automorphisms satisfying these conditions. Such 
automorphisms are called inner. Further call an inner automorphism unitary 
if c""'^ = c* (see (p9)) for the definition of the * operation). Wc can prove the 
following. 

Proposition 3. // T(x) = cxc~^ is a unitary inner automorphism and is 
representation of Hjy then there is a unitary operator Uc such that (j)(T(x)) = 
Uc4'{x)U~^ . Conversely for any unitary operator U on the representation space 
ofUN there is a unitary inner automorphism Tjj such that U(l){x)U~^ = 4>{Tij{x)) 
Thus there is a one-to-one correspondence between the set U(11n) of unitary 
inner automorphisms on C(Hjv) and quantum dynamics on the representative 
Hubert space. 
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This result is neither difficult nor surprising given the fact that the Hjv 
completely characterizes the kinematics of the system. It does however give us 
an alternative description and algebraic tools to study the dynamics. Thus we 
can study the effect of unitary operations on distribution functions[T2] using 
these transformations. Note however that we allow reducible representations 
now. The set of automorphisms of the group H^r is a subgroup of W(Hjv). 

In this work, we have concentrated on irreducible representations of H^r in 
which 7 acts maximally. By dropping the last assumption we can get all finite 
dimensional representations. The order of ^(7) in the representation (f> is the 
dimension. We can then use the products of these representations (actually we 
need some extra structure) to study unitary gates. We aim to explore this in 
future. Let us note some interesting relations in the case = 2". If u e H^r 
we will denote by (pk the representation in which 7^ =1. Let ai, i = 1,2,3 
denote the Pauli matrices and Ir the identity matrix of order r. Then 

02(^) = CcTiOCTi 02 ( ^) = crabs' 



where C 
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are a CNOT gate (C) and the phase gate {S) respectively [35]. We also note 
that in the general case 0„(X) is the cyclic shift operator. It can be efficiently 
constructed, for example, using full adder circuits with n + 1 ancillary qubits. 
Similarly, (pniZ) can be constructed using appropriate controlled phase gates as 
in the quantum Fourier transform. We also observe that iterating the simple 
relations (j>k{X'^) = 0fc_i(X) ^ I2 and 0fc(Z^) = I2 we obtain the 

interesting relations 

0„(X2') = 0„_fc(X) ® and MZ^") = h" ® ^n-k{Z) (57) 

These relations can be used to devise more efficient implementations. 

We conclude this section with a discussion of potential application of these 
constructions to quantum process tomography |36) . A quantum process is char- 
acterized by a completely positive map T acting on the operators on the system 
Hilbert space. If we have a complete set of phase-point operators {A{x, z)} then 
T is determined by its action on these. Let us assume that the dimension is odd 
so that we have a set of phase-point operators satisfying the full set of marginal 
conditions. Using Theorem [S]wc can prove the following. 

Proposition 4. Let T be quantum process (a CP map) given by 
T{A{x,z)) = Yj T{x',z' : x,z)A{x',z') 
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Here T{x',z' : x,z) is the "matrix" of T in the basis {A{x,z)} of phase-point 
operators. 



Then writing u = (1*1,^2) = {x,z) and u' = (u'j^jMj) = (x' , z') 

2 T(M'-u' : M-u) = (z' + '^\T (\z + ^><z + y l) 1^' + ^ 

(A/'-iu')i,(M-iu)i ^ ^ 

Here {\z}} and {\z'}} denote the ordered basis of eigenvectors of aM{Z) and 
aM'{Z) respectively. 

We do not prove it here since it is similar to the proof given in Theorem 
[5l Note that we are averaging over two indices now. To use the theorem we 
apply T to the projections \z}(z\ and measure the result in the basis H^')}- 
These transition probabilities yield the right hand side in the above equations. 
In principle, these equations can be inverted using the inverse Radon transforms 
(see section [5?T|) to yield the coefficients T{x', z' : x, z) and thus determining T 
(see [37| for a different perspective on phase-space tomography). Several opti- 
mizations are possible especially if we have some prior knowledge of the process. 
But we do not discuss these issues here as they merit a separate investigation. 

6 Discussion 

In this work we have analyzed quasiprobability distribution functions corre- 
sponding to quantum states. Our viewpoint is that these are the real coefficients 
of bases (or generally frames) in the space of hermitian operators. The choice of 
these bases is dictated by certain conditions we impose. This leads to expressing 
these bases or collection of phase-point operators in terms of operators represent- 
ing the Weyl-Heisenberg groups. In the language of frame theory [TUl [TT] these 
operators generate the Weyl-Heisenberg frames. We do not go into the intrica- 
cies of frame theory approach here. Our approach is more group-theoretic, em- 
phasizing the role of Weyl-Heisenberg groups in quantum kinematics. The other 
groups which play an important role are SL{2, Zjy) which yield the marginals. 
Conversely, we can use distribution functions to study these groups. We have 
given explicit formulas for the Radon transforms and their inversions. These 
can be used to solve the problem state or operator reconstruction. Even when 
we do not have sufficient data on marginals to invert the transforms we can 
get partial information about the state by taking generalized inverses |281 129j . 
We aim to address these and other issues on state and process estimation and 
reconstruction including the practical and computational aspects in future. 
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A Appendix 

Proof of Proposition [IJ We prove only (jlip . The invariance of scalar product 
and the standard measure on under rotation implies that 

jWcix.z: p)dx' = Jda;' J<e'("*+''^)>e-*("^+''^)dudi) 

= Jdx' J<e''("'*'+^^')>e-*("'=''+'''^'Mw'dw' 

= J<e'("'*'+^^')>e-'("'^'+^'"')(5(w')dw'dw' = J<e™'^'>e-'^'^'dw' 

Let \z'} be the eigenvectors of z' with eigenvalues z'. Then 

f dz'<e™''>= f dz'e'"'' (z'\p\z'} 

J a J a 

and equation pT|) follows. □ 
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Proof of Proposition [H We will prove the second formula only. The proof is 
similar for the first formula. Using the induced automorphism given in (|40p we 
get 

^ WoiM-'^x', M-^z' : p) = 



/V2 Zj Zj 



l-i5a,Ti + cn,o)(l-'56,T> + dn,o)((a+b)m+(c+d)n)^JV/4, ,(am + cn)(bm + dri)/2 

x'GZif m.n 



-San^ + a,^,o){l-Sb^^ + d^,o)((a+b)m + (c+d)nYN Ii^{am+cn)(bm+dn)l2 



4 s s-"- 

^-(a6m(m-l)/2 + cdn(ri-l)/2+(m + n)(cd+sgn(cd)/2))^^^^^j^-|m^^^^2')"^tj^"^'(5jj^g 
= J_^^;((c+d)")'Af/4^(cn)(dn)/2^-sgn(cd)ri/2-cd«(n-l)/2^^^^j-^^n^^-nz'^ 

Here we use that fact that cn, dn ¥= mod for any odd n since A^/w is even. 
Wc have to consider two cases separately suppose first that cd is odd. Then 

Af/2-1 



j r=0 



^(cd/2+j-^')(2r)_^ 



N/4-1 



2 (-l)4^^J+4^^Jc.M/2+J--')(2'-+i)+comp. conj. 



r=0 



= «z' - c<i/2|p|z' - cd/2} + (N/2 + z' ~ cd/2\p\N /2 + z' - cd/2))/2 

3 n odd 

ri<Ar/2 



We can prove the second case {cd odd) similarly. □ 
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